p 


AObl  206 


UNCLASSIFIED 


PENNSYLVANIA  STATE  UNIV  UNIVERSITY  PARK  APPLIED  RESE-- ETC  F/G  20/4 
A METHOD  FOR  CALCULATING  THE  FLOWFIELD  IN  THE  TAIL  REGION  OF  A — ETC(U) 
JUL  78  G H HOFFMAN  N00017-73-C-1418 

ARL/PSU/TM-78-211  NL 


DDC  FILE  CQPV  AD  AO  61  208 


A METHOD  FOR  CALCULATING  THE  FLOWFIELD  IN  THE  TAIL 
REGION  OF  A BODY  OF  REVOLUTION 

G.  H.  Hoffman 


Technical  Memorandum 
File  No.  TM  78-211 
July  19,  1978 

Contract  No.  N00017-73-C-1418 


Copy  No. 


L/h 


The  Pennsylvania  State  University 
Institute  for  Science  and  Engineering 
APPLIED  RESEARCH  LABORATORY^ 

Post  Office  Box  30 
— State  College,  Pa.  16801 


NAVY  DEPARTMENT 

NAVAL  SEA  SYSTEMS  COMMAND 

Approved  for  Public  Release 
Distribution  Unlimited 

7h 


o o if  o o 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dmta  Entered ) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3 PFCIPIFUT- 

17 


S CATALOG  NUMBER 


A METHOD  FOR  JIALCULATING  THE  JLOWFIELD  IN  THE 
TAIL  REGION  OF  A BPDY  OF  REVOLUTION. 


7.  AU  THORf a) 

)G.  H. /Hoffman 


9-  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Applied  Research  Laboratory 

P.  0.  Box  30 

State  College,  PA  16801 

It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Naval  Sea  Systems  Command 
Washington,  DC  20362 


s.  TYPE  OF  REPORT  & PERIOD  COVERED 

Technical  Memorandum 


6.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBERCs) 


J 1^317-73-0-1418? 


10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  & WORK  UNIT  NUMBERS 


f Ju 


14.  MONITORING  AGENCY  NAME  a ADDRESS (It  different  from  Controllint  Ottice)  15.  SECURITY  CLASS,  (ol  thla  report) 

’ / UNCLASSIFIED 


15a.  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 


16-  DISTRIBUTION  STATEMENT  (of  this  Report) 

Approved  for  Public  Release.  Distribution  Unlimited. 
Per  NAVSEA  - Oct.  16,  1978. 


[17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  If  different  from  Report) 


\ 18.  SUPPLEMENTARY  NOTES 


[ 19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  Identify  by  block  number) 


L. 


hydrodynamics 

flow 

boundary  layer 


ca  leu  I at  ion 


v streamline 

20.  ABSTRACT  (Continue  on  reverse  aide  If  necessary  and  Identify  by  block  number) 

(U)  A method  is  presented  for  calculating  a high  Reynolds  number  flow  field 
in  the  aft  region  of  a body  of  revolution  without  appendages  where  the 
boundary  layer  and  outer  inviscid  flow  strongly  interact.  The  method  is 
based  on  a correction  of  the  displacement  body  idea  to  effectively  represent 
streamline  curvature  effects  by  means  of  a simple  pressure  mapping.  The 
problem  is  solved  numerically  on  a digital  computer  by  iterating  for  the 
unknown  displacement  body  surface.  Typically,  from  6 to  8 iterations  arec^- 


DD  .ETn  1473  EOITION  OF  t NOV  65  1$  OBSOLETE 


unclassified 


91  OO'I 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Date  Entered 


SECURITY  CLASSIFICATION  OF  THIS  P AGEflWion  Data  Entered) 


Xrequired  for  convergence.  Available  computation  codes  are  used  for  the 
boundary  layer  flow  and  outer  inviscid  flow  while  a simple  integral  model 
is  used  for  the  wake.  Comparisons  of  computed  results  are  made  with 
experiment  and  other  numerical  solutions  for  three  bodies.  \ 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGErWTien  Data  Entered) 


Subject: 

References: 

Abstract: 


Acknow.T  edgement : 


A Method  for  Calculating  the  Flowfield  in  the  Tail  Region  of  a 
Body  of  Revolution 

See  Page  33 

A method  is  presented  for  calculating  a high  Reynolds  number  flow 
field  in  the  aft  region  of  a body  of  revolution  without  appendages 
where  the  boundary  layer  and  outer  inviscid  flow  strongly  interact. 
The  method  is  based  on  a correction  of  the  displacement  body  idea 
to  effectively  represent  streamline  curvature  effects  by  means  of 
a simple  pressure  mapping.  The  problem  is  solved  numerically  on 
a digital  computer  by  iterating  for  the  unknown  displacement  body 
surface.  Typically,  from  6 to  8 iterations  are  required  for 
convergence.  Available  computation  codes  are  used  for  the  boundary 
layer  flow  and  outer  inviscid  flow  while  a simple  integral  model 
is  used  for  the  wake.  Comparisons  of  computed  results  are  made 
with  experiment  and  other  numerical  solutions  for  three  bodies. 

This  work  was  sponsored  by  Naval  Sea  Systems  Command,  Codes  NSEA-0351 
and  NSEA-03133 . 


ACCESSION  for 

, 

NT  IS 

W'ri'e  Section  frf 

DDC 

B.ff  Section  □ 

UNANNOUNCED 

□ 

JUST  r ICAT 

ION 

— 

f BY 

OOTJMWABril 

1 CKtS 

Out 

-2-  July  19,  1978 

GHH: jrp 

Table  of  Contents  page 

Abstract  1 

Acknowledgement  1 

Nomenclature  3 

List  of  Figures 5 

List  of  Tables  7 

Introduction  8 

Pressure  Mapping  Hypothesis  10 

Strong  Interaction  Calculation  Method  14 

Wake  Model  17 

Calculated  Results  21 

Lyon  Model  A Body 23 

Modified  Spheroid  25 

F-57  Body  27 

Conclusions 30 

References  33 

Figures 35 


Nomenclature 


reference  area 


o 

drag  coefficient  (=  drag/'ipU^  A^) 


2 

wall  friction  coefficient  (=  x /*spu  ) 


pressure  coefficient 


axisymmetric  shape  factor 


body  length 


distance  normal  to  stream  surface 


static  pressure 


velocity  along  streamline 


radial  distance 


local  body  radius 


radius  of  curvature 


Reynolds  number  (=  U L/v) 


distance  along  stream  surface 


velocity  in  x-direction 


inviscid  velocity  at  r = r. 


free  stream  velocity 


distance  along  body  surface 


axial  distance 


distance  normal  to  body  surface 


boundary-layer  thickness 
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axisymmetric  physical  displacement  thickness 


mass  deficit  area  (=  \ (1  — — ) rdy) 

J 


eddy  viscosity 


V 


kinematic  viscosity 
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Nomenclature  (Cont.) 


<t> 

P 

0 

T 


W 


angle  between  tangent  to  body  surface  and  x^-direction 

density  „ 

f u 

momentum  deficit  area  (=  \ — (1-  — ) rdy) 

J u u 

0 e e 

surface  shear  stress 


All  other  quantities  and  subscripts  are  defined  in  the  text. 

All  quantities  in  the  text  (unless  otherwise  specified)  are  made  dimension- 
less as  follows: 

distance  with  respect  to  L 

2 

area  with  respect  to  L 
velocity  with  respect  to  U 

2 

pressure  with  respect  to  pU 
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Introduction 

The  most  difficult  part  of  the  flow  field  to  compute  about  a body  of 
revolution  is  in  the  vicinity  of  the  tail-wake  juncture,  referred  to  as 
the  strong  interaction  region.  For  applications  this  region  is  important 
in  two  respects.  First,  the  drag  depends  to  a large  extent  on  the  flow  in 
this  region  and  second,  when  a propeller  or  other  propulsive  device  is  added 
it  always  operates  in  the  thick  boundary  layer  near  the  tail. 

This  report  addresses  the  problem  of  calculating  the  flow  field  in  the 
strong  interaction  region  of  a body  of  revolution  at  zero  angle  of  attack 
where  the  boundary  layer  does  not  separate.  The  aim  has  been  to  develop  an 
engineering  tool  capable  of  predicting  this  flow  field  with  reasonable 
accuracy  using  a blend  of  available  calculation  methods  plus  a new  twist  to 
the  displacement  body  idea.  Only  the  incompressible  case  is  considered. 

The  phenomenon  that  occurs  in  the  tail-wake  juncture  region  of  a body 
of  revolution  can  be  explained  as  follows:  Toward  the  aft  end  of  a body 
with  a pointed  tail  the  boundary-layer  thickness  becomes  comparable  to  the 
body  radius  thereby  significantly  displacing  the  streamlines  of  the  outer 
inviscid  flow  which  in  turn  modifies  the  body  pressure  distribution.  In 
addition,  the  streamlines  of  the  outer  flow,  sensing  the  wake  line  of 
symmetry  approaching,  undergo  a significant  change  in  direction  over  a 
relatively  short  distance  thereby  impressing  a non-negligible  normal 
pressure  gradient  on  the  boundary  layer.  To  treat  the  resulting  strong 
interaction,  the  equations  governing  the  boundary  layer  and  the  outer 
flow  need  to  be  solved  simultaneously  and  the  assumptions  of  conventional 
boundary-layer  theory  modified. 


Over  the  past  10  to  15  years  a number  of  methods  have  appeared  for 
solving  the  strong  interaction  problem.  The  most  common  method  is  the 
displacement  body  approach,  the  category  under  which  the  present  work 
falls.  In  this  method  the  boundary-layer  displacement  thickness  is  added 
to  the  original  bodv  to  represent  the  thickening  effect  that  the  outer 
flow  senses  and  the  boundary  layer  calculated  as  though  no  normal  pressure 
gradient  existed.  The  new  aspect  of  the  present  work  is  that  the  displacement 
body  idea  is  modified  by  a simple  pressure  mapping  which  effectively  represents 
streamline  curvature  effects  on  the  boundary  layer. 

A comprehensive  review  up  to  1976  on  methods  of  treating  the  strong 
interaction  problem  has  been  given  by  Huang  et.  al.  [1].  Since  that  time 
two  other  papers  dealing  with  new  approaches  to  the  problem  have  appeared. 

The  work  by  Schetz  and  Favin  [2]  is  the  most  exact  approach  to  the  problem 
to  date.  They  numerically  solve  the  Navier-Stokes  equations  with  a one 
equation  transport-type  turbulence  model  in  a truncated  region  which  includes 
the  aft  portion  of  the  body  and  part  of  the  wake.  The  effect  of  a propeller 
positioned  at  the  tail  is  modeled  by  an  actuator  disk.  Unfortunately  no 
comparison  with  experiment  is  given.  The  recent  paper  by  Dyne  [3]  is  a 
streamline  curvature  method  in  which  the  effects  of  turbulence  in  the 
boundary  layer  have  been  accounted  for  approximately.  The  calculation 
region  where  the  streamline  curvature  method  is  applied  is  the  aft  portion 
of  the  body  and  the  near  wake,  similar  to  the  computational  region  of 
Schetz  and  Favin.  The  virtue  of  the  streamline  curvature  method  is  that 
the  normal  pressure  gradient  caused  by  streamline  curvature  is  taken  into 
account  while  the  simplicity  of  a boundary  layer-like  formulation  is  re- 
tained. Agreement  with  experiment  for  the  one  case  presented,  a modified 
spheroid,  is  good. 
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2.  Pressure  Mapping  Hypothesis 

In  the  usual  displacement  body  treatment,  the  pressure  from  that 
body  is  referred  back  to  the  original  body  along  a boundary-layer  normal 
coordinate.  As  Weinbaum  et  al.  [4]  point  out,  this  transfer  implies  that 
the  pressure  gradients  on  the  two  bodies  are  related  by 
d£  = R*  d£ 

ds  R ds*.  u; 

w w 

where  the  subscript  w denotes  the  surface  of  the  original  body,  the 
superscript  * the  surface  of  the  displacement  body,  s is  arc  length  and 
R is  the  local  radius  of  curvature.  Equation  (1)  is  a consequence  of 
the  differential  arc  length  relation  between  the  displacement  body  and 
the  original  body, 

~ dsw  = ds*,  (2) 

w 

which  holds  provided  the  bodies  are  sufficiently  smooth  and  the  curvatures 
are  changing  slowly. 

The  two-dimensional  analysis  of  Weinbaum  et  al.  for  a thick  viscous 
layer  at  intermediate  Reynolds  numbers  will  now  be  extended  to  the  axisym- 
metric  case  at  high  Reynolds  numbers.  The  object  of  this  analysis  is  to 
determine  how  the  conventional  displacement  body  idea  must  be  modified  to 
account  for  centrifugal  effects  in  a turbulent  wall  shear  layer  at  high 
Reynolds  numbers.  With  the  assumption  that  the  turbulence  field  can  be 
modeled  by  an  eddy  viscosity,  the  appropriate  set  of  equations  in  mean 


flow  intrinsic  coordinates  is, 


9 If"  = " t(1+e>r 


= iE. 

R 3n  * 


(4) 
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These  equations  are  the  axisymmetric  turbulent  analogs  of  Eqs.  (2)  and 
(3)  of  Weinbaum  et  al.  The  above  equations  have  been  made  dimensionless 


with  velocities  referred  to  free  stream  velocity  U , pressure  to 
2 


pU^  , lengths  to  reference  length  L and  eddy  viscosity  to  kinematic 
viscosity  V.  Equations  (3)  and  (4)  are  derived  by  considering  the  full 
axisymmetric  Navier-Stokes  equations  in  general  orthogonal  coordinates, 
introducing  the  usual  Reynolds  turbulence  -averaging,  making  the  boundary- 
layer  approximation  and  finally,  specializing  the  result  to  intrinsic 
coordinates  of  the  mean  turbulent  flow. 

Equations  (3)  and  (4)  are  now  ordered  in  terms  of  the  Reynolds 
number  by  stretching  the  normal  coordinate  according  to  the  usual  high 
Reynolds  number  boundary-layer  scaling: 


Re 


* 


(5) 


In  the  intermediate  Reynolds  number  case  of  Weinbaum  et  al.  both  spatial 
variables  are  scaled  by  the  inverse  square  root  of  the  Reynolds  number. 
With  the  scaling  given  by  Eq.  (5)  the  equations  of  the  thick  shear  layer 
become 


If  + f lsm«)rU]. 


3s  r 3n 


3tv 


2 

= 

R 


Re'2  & 
Ke  3n 


(6) 

(7) 


The  convective  acceleration  term  is  eliminated  by  differentiating  Eq.  (7) 
with  respect  to  s,  then  substituting  the  result  into  Eq.  (6)  which  gives 
the  following  equation: 


iL  ,iE)  + A iE  = _ _L_  ‘ 

9"  (9s>  Re^R  95  Re*5  <R> 


3R 


3s 


Re  2Rr 


I?  i<i+e>r  Hi- 


(8) 
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Note  that  Eq.  (8)  is  the  high  Reynolds  number  axisymmetric  analog  of 
Eq.  (5)  of  Weinbaum  et  al.  and  describes  the  variation  of  streamwise 
pressure  gradient  across  the  shear  layer. 

To  simplify  the  integration  of  the  left-hand  side  of  Eq.  (8)  the 
local  radius  of  curvature  is  approximated  by 

R = R + — V . 

W Re55 


(9) 


Then  integration  of  Eq.  (8)  from  the  wall  (n  = 0)  to  a point  in  the  shear 
layer  gives  the  streamwise  pressure  gradient  at  n as 


3£  = ,Rw  8£ 

3s  ''R  ’ ^3 s J 


n n 

1 (1  ,dR.  [ 2 - 2 [ R 8 , 

- 75  {7  <dl>  ’ <*”  - ~2  V 3B 

w Re  R w J R J 


3=[(l+e)r^]  dn}. 
3n 


Equation  (10)  shows  that  the  change  in  streamwise  pressure  gradient  normal 
to  the  wall  is  dominated  by  streamline  curvature  effects,  represented  by 

-h 

the  first  term.  The  remaining  terms,  0(Re  ),  represent  respectively 
changes  due  to  streamwise  variation  in  wall  curvature  and  laminar/turbulent 
diffusion. 

One  next  applies  Eq.  (10)  at  the  surface  of  the  displacement 
body,  assuming  that  body  to  be  a streamline  of  the  actual  flow  (which  in 
reality  it  is  not) . Thus  the  streamwise  pressure  gradient  along  the 
displacement  body  is 


l£  = 

3s* 


R 2 

(— ) 
VR*' 


^3s' 


0(Re-i5). 


w 


(11) 


The  pressure  gradient  on  the  displacement  body  is  now  transferred  back  to  the 
original  body  along  the  boundary-layer  normal  according  to  Eq.  (1)  which 
gives  the  result 


wn 


<{r>  (H>  + 0(Re"‘i)- 

w w 


(12) 


(10) 
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where  the  subscript  wn  denotes  the  wall  pressure  gradient  In  the 
presence  of  a normal  pressure  field.  If  the  following  arc  length  stretch- 
ing is  now  defined: 

ds  = (^— ) ds  — = ds*,  (13) 

w R wn 
w 

then  Eq.  (12)  becomes 

|f.  -(|f)  <14) 

w 

which  is  the  result  obtained  by  Weinbaum  et  al.  Equations  (13)  and  (14) 
are  the  essence  of  Weinbaum' s hypothesis  which  states:  "The  pressure 
gradient  along  the  displacement  body  constructed  in  the  conventional 
manner  using  first-order  boundary-layer  theory  will  be  the  same  as  in  a 
viscous/ turbulent  shear  layer  theory  which  includes  centrifugal  effects 
provided  the  pressure  on  the  displacement  body  is  referred  back  to  the 
original  body  in  a manner  that  preserves  the  arc  length."  The  present 
analysis  shows  that  Weinbaum' s hypothesis  is  unchanged  by  either  axial 
symmetry  or  a turbulent  eddy  viscosity. 

In  regions  of  large  changes  in  the  streamline  curvature  both  neglected  terms  in 
Eq.  (10)  may  become  0(1)  which  will  invalidate  Weinbaum' s hypothesis.  In  the  high 
Reynolds  number  case  the  diffusion  term  becomes  equal  in  importance  to  the  term 
arising  from  streamwise  variations  in  the  wall  curvature,  in  contrast  to  the  inter- 
mediate Reynolds  number  case  where  the  diffusion  term  is  0(Re  ^) . This  behavior  is 
a consequence  of  the  different  scalings. 

The  crucial  assumption  in  the  preceding  development  and  the  one  directly  respon- 
sible for  the  simple  form  of  Weinbaum' s hypothesis  is  Eq.  (9),  the  approximation  of 
the  local  radius  of  curvature.  Physically  this  approximation  means  that  the  stream- 
lines are  locally  concentric  which  is  usually  not  the  case  near  the  tail  of  the  body. 
The  validity  of  Weinbaum' s hypothesis  with  its  simplifying  assumptions  can  only  be 
determined  by  comparison  with  experiment  over  a range  of  strong  interaction  cases. 
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3.  Strong  Interaction  Calculation  Method 

The  solution  for  the  location  of  the  displacement  body  requires 
i that  the  inviscid  pressure  distribution  be  known  at  the  surface  of  the 

original  body.  In  turn,  the  solution  for  the  pressure  distribution  on 
the  original  body  requires  that  the  displacement  body  be  known.  Thus  a 
coupled  nonlinear  boundary  value  problem  exists  which  must  be  solved  by 
iteration. 

The  iteration  procedure  used  here  follows  closely  the  schemes  used 
in  displacement  body  calculations  found  in  the  literature  [1,  5,  6].  The 
general  steps  in  the  iteration  are  as  follows: 

First  iteration 

1.  Compute  the  inviscid  pressure  distribution  for  the  basic  body. 

2.  Modify  the  pressure  distribution  from  step  1 on  the  afterbody  to 
eliminate  the  rear  stagnation  point,  it  is  exists,  so  that  boundary- 
layer  separation  will  be  avoided.  Guess  the  wake  pressure  distribution 
for  1 < Xq  < 2.  Step  2 gives  the  first  approximation  of  the  pressure 
distribution  over  the  body-wake  combination. 

3.  With  the  pressure  distribution  from  step  2 compute  the  boundary 
layer  along  the  actual  body  and  the  wake  at  least  one  body  length 
downstream  of  the  tail.  This  step  gives  the  first  approximation 
to  the  displacement  body. 

Second  and  succeeding  iterations 

1.  Compute  the  inviscid  pressure  distribution  using  the  displacement 
body  from  the  last  step  of  the  previous  iteration. 

2.  Transfer  the  pressure  distribution  from  the  displacement  body  back  to 


the  basic  body  and  the  wake  centerline,  preserving  arc  length  (Weinbaum's 
hypothesis)  . 
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3.  With  the  pressure  distribution  from  step  2,  compute  the  boundary 

layer  along  the  actual  body  and  the  wake  to  obtain  a new  estimate  of  the 
displacement  body  shape. 

The  iteration  cycle  is  continued  until  convergence  is  achieved.  The 
pressure,  being  the  most  sensitive  variable,  is  monitored  until  its 
maximum  change  over  successive  iterations  is  sufficiently  small.  The 
steps  as  outlined  above  will  be  elaborated  upon  as  each  part  of  the 
calculation  cycle  is  taken  up  in  detail.  What  distinguishes  the  present 
displacement  body  method  from  others  is  the  use  of  Weinbaum's  hypothesis 
to  account  for  centrifugal  effects. 

The  philosophy  followed  here  has  been  to  use  available  calculation 
methods  as  much  as  possible  for  the  components  of  the  strong  interaction 
calculation  - the  inviscid  flow,  the  boundary  layer  and  the  wake.  For  the 
inviscid  flow  the  method  of  Hess  and  Smith  [7]  is  used.  For  the  boundary- 
layer  calculation  transverse  curvature  terms  are  included  in  the  equations 
(the  so-called  thick  boundary-layer  approximation).  These  equations  are 
solved  by  a modified  Cebici-Smith  finite  difference  approach  [8]  which 
makes  use  of  the  second-order  accurate  form  of  Keller's  box  method.  The 
code  in  use  was  generated  at  ARL  and  includes  a two-layer  algebraic 
turbulence  model  composed  of  the  Cebici  form  for  the  outer  eddy  viscosity 
[8]  and  the  Crawf ord-Kays  form  for  the  inner  eddy  viscosity  [9].  The  inner 
eddy  viscosity  has  been  modified  for  a thick  axisymmetric  boundary  layer 
and  reformulated  to  work  properly  in  the  region  of  a pointed  tail.  For 
the  wake  calculation  an  integral  method  is  used  which  is  covered  in  detail 
in  the  next  section. 

In  the  axisymmetric  boundary  layer  three  definitions  of  the  displacement 
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thickness  exist.  The  one  appropriate  to  the  displacement  body  is  the 
so-called  physical  definition  S^*  which  is  derived  from  the  mass  flow 
defect  in  the  annular  boundary  layer  wrapped  around  the  body: 


rn+S.,  *cos<f>  rQ+6cos<f> 


iO,v,l  >0 

\ 2nr  ugdr  = I 


2iTr  (u  -u)dr. 
e 


Integration  of  the  above  equation  with  the  aid  of  the  geometric  relation 


r = Tq  + ycos<{>  in  the  right-hand  integral  yields 


0 

1 

0 


S^*r  + ^i*^cos<J>  ~ l (1  - )rdy  = A*, 

1 ue 


(15) 


where  A*  is  the  mass  deficit  area.  From  Eq.  (15)  the  solution  for  6^* 
is 


h 


61*  = 


(2A *costp  + rQ  ) - rQ 

C0S(j) 


(16) 


From  the  geometry  shown  in  Fig.  1 the  coordinates  of  the  displacement  body 
are  given  by: 


(r0}6*  = r0  + 6i*cos<*>’ 

^s*  = xo  “ 6i*sin<J,» 


or  finally, 

2 ** 

(Vs*  = (2A*cos<^  + r0  ) « (17) 

2 ^ 

(xq)s*  = xo  “ tan<J)[(2A*cos<|)  + rQ  ) - rQ] . 


(18) 
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4 . Wake  Model 

Because  of  the  presence  of  a wake  the  displacement  body  continues 

to  infinity  downstream  of  the  end  of  the  actual  body,  but  its  importance 

aft  of  the  tail  in  establishing  curvature  effects  on  the  body  pressure 

distribution  is  mainly  in  the  near  wake  region  (1.0  < < 1.1).  Thus 

the  wake  portion  of  the  displacement  body  must  be  included  in  a strong 

interaction  calculation  but  the  calculation  of  the  wake  can  be  terminated 

about  one  body  length  downstream  of  the  tail  without  producing  any  adverse 

effects  in  the  overall  solution. 

The  wake  is  modeled  in  the  simplest  manner  possible  by  using  an 

integral  approach.  The  two  variations  in  this  type  of  treatment  involve 

either  using  an  entrainment  relation,  as  in  the  work  of  Nakayama,  Patel 

and  Landweber  [10],  or  using  a relation  between  shape  factor  and  inviscid 

velocity,  as  in  the  work  of  Granville  [11].  Granville's  variant  was  chosen 

because  it  leads  to  a simple  closed  form  solution  whereas  the  entrainment 

method  must  be  integrated  numerically.  Furthermore,  as  Patel  and  Guven 

point  out  [12],  the  assumed  variation  of  H with  u is  not  crucial  to 

e 

the  determination  of  the  asymptotic  wake  behavior.  But,  as  will  be  seen, 
the  asymptotic  behavior  does  strongly  depend  on  the  initial  values  of  u 

e 

and  0. 


In  the  wake  the  axisymmetric  integral  momentum  equation  with  normal 


pressure  variations  neglected  reduces  to 
...  d£nu 

£ + <H+2> 
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By  curve  fitting  axisymmetric  wake  data  Granville  finds  the  following 
relation  between  H and  u [11]: 

6 1/q 

£nu 

H = 1 + (H  - 1)  - ) (21) 

t u„ 

et 

where  the  subscript  t denotes  the  tail  of  the  body,  = 1,  and  q 

is  found  by  Granville  to  be  7.  With  the  aid  of  Eq . (21),  Eq.  (20)  can 
be  integrated  giving  the  result 


6 


where 


0 = 0 u 

co  e 


n 2q  + 3 + qH 

3 " ' 1 T q 


(22) 


(23) 


and  0^  is  the  momentum  deficit  area  infinitely  far  downstream  in  the 
wake.  In  the  present  context  ug  in  Eq.  (22)  is  the  inviscid  velocity 
on  the  wake  centerline  obtained  by  transferring  the  pressure  on  the  dis- 
placement body  to  the  wake  centerline  using  Weinbaum's  hypothesis. 

Since  we  assume  that  Eq.  (22)  holds  at  the  initial  wake  station,  we 
can  determine  0^  in  terms  of  0^,  known  from  the  boundary-layer  calculation 
at  the  tail,  and  ug  , known  from  the  pressure  mapping.  Thus 


0_ 


0^uf 


t et  ’ 

and  hence  at  any  other  station 
H by  Eq.  (21).  Note  that  Eq. 
drag  formula  is  derived. 

With  0 and  H known  at 
area  A*  is  given  by 


(24) 

in  the  wake  0 is  given  by  Eq.  (22)  and 
(24)  is  the  expression  from  which  Granville's 

a particular  wake  station,  the  displacement 


A*  = H0, 


(25) 
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and  finally,  the  radius  of  the  displacement  body  is  given  by  Eq.  (17) 
which  in  the  wake  reduces  to 

(r  ) = (2A*)*5.  (26) 

u 6* 

The  initial  values  of  the  wake  parameters  to  be  used  in  Eq.  (24)  are 
obtained  by  requiring  r^  and  H to  be  continuous  at  the  body-wake 
juncture.  Obviously  the  displacement  body  radius  must  be  continuous  for 
a smooth  pressure  distribution  to  exist.  The  requirement  of  continuous 
H is  based  on  the  observation  that  a non-zero  tail  angle,  which  causes 
a mismatch  or  overlap  in  boundary-layer  coordinates  with  wake  (cylindrical) 
coordinates,  does  not  appear  to  introduce  a discontinuity  in  the  measured 
shape  factor.  These  two  conditions  are  sufficient  to  determine  the  remaining 
parameters  at  the  body-wake  juncture.  With  r^  and  H known  at  the 
initial  wake  station  t,  and  with  uet  given  by  the  pressure  mapping, 
the  other  parameters  are  given  by: 

A\  = i (r6*^  ’ (27) 

and 

Q,  = A*t  (28) 

Ht 

Although  rg*  and  H are  forced  to  be  continuous  at  = 1.0,  the 
other  parameters  at  that  point  may  be  discontinuous  since  a finite 
difference  boundary-layer  solution  is  being  joined  to  an  integral  wake 
solution.  These  discontinuities  have  been  found  to  be  quite  small.  Other 
ways  of  treating  the  boundary  layer-wake  juncture  were  tried  but  none  gave 
as  satisfactory  results  as  the  simple  treatment  just  described. 

The  Mangier  transformation  used  in  the  Cebici-Smith  boundary-layer 
method  introduces  a singularity  into  the  equations  of  motion  at  the  tail 
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when  the  body  radius  there  is  zero.  Consequently  the  boundary  layer 
solution  cannot  be  calculated  all  the  way  to  the  tail  with  this  method. 
This  difficulty  is  overcome  by  computing  the  boundary-layer  solution  to 
about  Xq  = 0.99,  then  obtaining  values  of  r^  and  H at  the  tail 
station  by  parabolic  extrapolation. 


L 
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5 . Calculated  Results 

To  begin  a strong  interaction  calculation,  the  body  pressure 
distribution  obtained  by  solving  the  potential  flow  problem  about  the  basic 
body  without  a displacement  thickness  added  is  modified  in  the  tail  region 
to  approximate  the  converged  strong  interaction  solution.  The  modified 
pressure  distribution  is  then  extrapolated  into  the  wake.  The  modification 
of  the  basic  pressure  distribution  begins  at  about  Xq  = 0.85  with  the 
modified  curve  passing  smoothly  through  Cp  = 0.2  at  Xq  = 1,0,  then 
decaying  rapidly  to  zero  in  the  wake.  The  wake  Cp  distribution  is 
approximated  by 

Cp  = 0.2  sech  [ 20. 6<Xq-1) ] , xQ  >.  1.  (29) 

Even  though  the  initial  guess  is  a reasonably  good  one  the  second 
iteration  usually  produces  an  overcorrection  to  the  pressure  distribution 
which  if  carried  further  leads  to  divergence.  This  type  of  behavior  was 
also  found  by  Nakayama  et  al.  [10].  Their  solution  to  the  overcorrection 
problem  was  to  compute  a pressure  distribution  for  the  second  iteration 
by  averaging  the  initial  guess  (first  iteration)  and  the  pressure  distribution 
obtained  from  the  resulting  first  iteration  boundary-layer  solution.  After 
the  second  iteration,  averaging  was  found  to  be  unnecessary  and  a convergent 
procedure  followed.  This  same  procedure  was  found  to  work  successfully  in 
the  present  calculation  method.  The  only  difference  is  that  here  the 
averaging  is  performed  after  the  application  of  the  pressure  mapping. 

During  the  initial  stages  of  this  work  convergence  was  found  to  be 
slowest  in  the  region  near  the  tail-wake  juncture  because  spurious  jumps 
in  the  pressure  distribution  occurred.  These  jumps  were  found  to  be 
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caused  by  discontinuities  in  the  slope  and  curvature  of  the  displacement 
body  at  the  juncture  point  which  in  turn  stemmed  from  the  incompatibility 
of  a finite  difference  boundary  layer  solution  patched  to  an  integral  wake 
solution. 

Huang  et  al.  [1],  whose  strong  interaction  calculation  method  the 
present  one  resembles,  also  encountered  spurious  pressure  jumps  in  the 
tail-wake  juncture  region  where  they  represented  the  displacement  body  by 
a polynomial  between  0.95  < x^  < 1.05.  When  a fifth-order  polynomial  was 
used  to  ensure  continuity  of  radius,  slope  and  curvature  at  the  juncture 
points  the  pressure  jumps  disappeared.  Their  polynomial  fairing  thus 
provided  a smoothing  of  the  pressure  in  the  tail-wake  juncture  region. 

Smoothing  in  the  present  calculations  was  performed  when  needed  by 
fairing  out  the  spurious  jumps  in  the  pressure  distribution  by  hand  rather 
than  fairing  the  displacement  body  tail-wake  juncture.  Although  the 
present  method  of  eliminating  the  pressure  jumps  is  primitive  and  involves 
judgement  on  the  part  of  the  user,  the  iteration  process  does  converge, 
usually  within  6 to  8 iterations,  and  the  pressure  jumps  disappear  by  the 
final  iteration.  Convergence  is  considered  achieved  when  the  maximum 
change  in  u^  (rather  than  C^)  is  less  than  about  one  percent.  The 
maximum  change  always  occurs  close  to  the  tail-wake  juncture. 

Present  experience  with  the  Cebici-Smith  boundary-layer  calculation 
method  shows  that  even  though  the  sublayer-buffer  region  of  the  turbulent 
boundary  layer  may  be  adequately  resolved  at,  say,  x^  = 0.95,  it  will  not 
be  at  Xq  = 0.99.  The  fault  lies  with  the  Mangier  transformation  which  is 
singular  at  the  tail  when  body  radius  is  zero.  The  step  size  problem  may 
be  demonstrated  as  follows:  For  the  first  step  normal  to  the  body  in  the 
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boundary  layer  the  physical  step  size  Ay^  Is  related  to  the  transformed 
step  size  An^  by  (see  Ref.  8,  p.  260) 

Ani 

Ay,  = constant  * 

1 r 


so  that  as  r^  -*•  0 with  Ari^  fixed,  the  physical  step  size  in  the 
boundary  layer  becomes  large  even  though  Ari^  may  be  extremely  small. 

Solutions  were  computed  by  the  present  method  for  three  axisymmetric 
bodies  for  which  experimental  data  have  been  published.  The  three  cases 
are  the  model  A body  of  Lyon  [13],  a modified  spheroid  [14]  and  a low- 
drag  body  [15],  The  first  two  bodies  have  non-zero  tail  angles  while  the 
last  is  nearly  cusped.  In  each  case  comparisons  with  the  data  are  given 
for  body  pressure  coefficient  C^,  momentum  deficit  area  0,  shape  factor  H 
skin  friction  coefficient  plus  selected  velocity  profiles.  Where 

possible  comparisons  are  also  made  with  published  theoretical  results. 


and 


Lyon  Model  A Body 

Calculations  were  made  at  a body  chord  Reynolds  number  of  2.09  x 10^, 
the  same  as  in  Lyons'  experiments  [13].  The  transition  location  in  the 
calculation  was  adjusted  so  that  0 in  the  laminar-turbulent  region 
(0.1  < Xq  < 0.4)  closely  fitted  the  experimental  data  given  in  Fig.  24 
of  ref.  13.  These  data  correspond  to  the  case  with  a screen  ahead  of  the 
model.  The  turbulent  eddy  viscosity  in  the  calculation  is  activated 
gradually  downstream  of  the  transition  station  by  using  the  intermittency 
factor  of  Chen  and  Thyson  [16]. 

The  geometry  of  the  Lyon  model  A body  is  shown  in  Fig.  2.  Figure  3 
presents  a comparison  between  the  calculated  and  measured  pressure  distrib- 
utions on  the  aft  portion  of  the  body.  The  calculated  pressure  distribution 
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along  the  wake  centerline  is  also  shown  where  no  measurements  were  made. 

In  addition,  the  potential  solution  (no  boundary-layer  interaction)  and  the 
strong  interaction  solution  of  Wakayama,  Patel  and  Landweber  [10]  are 
presented  for  comparison.  The  agreement  of  the  solution  of  Wakayama  et  al. 
with  the  data  is  slightly  better  than  the  present  solution.  In  Fig.  4 
the  momentum  deficit  area  distribution  as  predicted  by  the  present  method 
and  that  of  Nakayama  et  al.  is  compared  with  the  experimental  data.  On  the 
body  both  predictions  agree  closely  with  experiment.  In  the  wake  the 
calculated  values  of  Nakayama  et  al.  are  above  those  of  the  present  method. 

Since  a strong  interaction  calculation  necessarily  includes  the  wake, 
one  of  the  results  is  the  body  drag  which  is  related  to  the  momentum  area 
at  downstream  infinity  by 


where  is  the  reference  area.  A comparison  between  experiment  and  several 

calculations  of  CQ  (based  on  body  surface  area)  for  the  Lyon  model  A body 
is  given  in  Table  1. 


Source 

CD 

Experiment  - Lyon  [13] 

.00436 

Present  calculation 

.00401 

Calculation  - Nakayama  et  al.  [10] 

.00410 

Calculation  - Myring  [5] 

.00397 

Table  1.  Drag  Coefficient  for  Lyon  Model  A Body. 

The  spread  in  the  three  calculated  values  is  only  3.2  percent. 

The  results  for  and  H are  shown  in  Figs.  5 and  6 respectively. 

Lyons'  measured  values  of  C^.  display  an  erroneous  trend  near  the  tail  of 
the  body  but  agree  well  with  the  present  calculation  over  the  front  portion. 


1 
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The  shape  factor  Is  predicted  better  over  the  last  30  percent  of  the  body 
by  the  integral  method  of  Nakayama  et  al.  than  ti.e  present  finite  difference 
method.  The  reason  for  the  poor  showing  of  the  present  method  for  H is 
attributed  to  inadequency  of  the  algebraic  turbulence  model  in  the  tail 
region  even  though  transverse  curvature  effects  are  included  in  the  model 
as  well  as  in  the  boundary-layer  equations.  Velocity  profiles  are  shown 
at  Xg  = 0.6  in  Fig.  7,  where  the  thin  boundary-layer  approximation  is 
valid,  and  at  x^  = 0.9  in  Fig.  8 in  the  strong  interaction  region  where 
transverse  curvature  effects  are  important.  At  x^  = 0.6  agreement  with 
experiment  is  excellent  but  at  xq  = 0.9  the  shape  of  the  predicted  profile 
is  too  full  near  the  surface.  The  latter  behavior  is  typical  of  the 
present  method  for  the  three  cases  presented  herein. 

Modified  Spheroid 

The  body  used  by  Patel,  Nakayama  and  Damian  [14]  in  their  experiments 
was  a spheroid  of  approximately  6:1  fineness  ratio  with  the  rear  7 percent 
removed  and  a tangent  cone  added  to  prevent  boundary-layer  separation.  The 
modified  spheroid  geometry  is  shown  in  Fig.  9.  Corresponding  to  the 
experiment  of  Patel  et  al.,  strong  interaction  calculations  have  been 
made  at  a chord  Reynolds  number  of  1.262  x 10^  and  the  eddy  viscosity 
switched  on  at  x^  = 0.05  to  simulate  the  boundary-layer  trip  on  the  actual 
body.  Comparisons  between  present  calculations,  those  of  Nakayama  et  al. 

[10]  and  the  experiments  of  Patel  et  al.  are  given  in  Figs.  10  - 15.  The 
pressure  distribution  on  the  body  and  on  the  wake  centerline  as  well  as  the 
momentum  deficit  area  distribution,  shown  in  Figs.  10  and  11  respectively, 
are  predicted  well  by  both  methods  except  that  the  calculated  pressure  peak 
near  the  tail  is  slightly  high.  Note  that  for  this  body  the  departure 
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from  the  potential  flow  pressure  distribution  begins  at  xQ  = 0.70.  For 
both  sets  of  calculations  the  skin  friction  coefficient  is  generally  higher 
than  the  experimental  values  as  seen  in  Fig.  12.  The  present  calculation 
indicates  that  the  boundary-layer  nearly  separates  at  x^  = 0.94  but  then  recovers. 
This  trend  is  not  shown  by  the  integral  method  result  of  Nakayama  et  al. 

The  boundary-layer  shape  factor  prediction  by  the  present 
method,  shown  in  Fig.  13,  becomes  increasingly  less  accurate  as  the  tail  is 
approached,  as  occurred  for  the  Lyon  body.  The  integral  method  calculation 
of  Nakayama  et  al,  is  closer  to  the  experimental  values  in  the  vicinity 
of  the  tail  but  in  the  wake  the  present  method  gives  a better  prediction 
of  H.  Nakayaiaa  et  al.  point  out  concerning  their  calculations  for  the 
modified  spheroid  that  lower  values  of  H and  higher  values  of  have 

a mutually  cancelling  effect  in  the  momentum  equation  in  an  adverse  pressure 
gradient  which  explains  why  their  method  as  well  as  the  present  one  predict 
the  momentum  deficit  area  well  near  the  tail. 

The  predicted  velocity  profiles  for  the  modified  spheroid,  shown  in 
Figs.  14  and  15,  behave  in  much  the  same  way  as  for  the  Lyon  body.  At 
Xq  = 0.662  the  prediction  agrees  well  with  experiment  but  at  x^  = 0.96, 
in  the  strong  interaction  region,  the  predicted  velocities  are  too  high 
near  the  wall. 

Finally,  the  experimental  and  calculated  drag  coefficients  for  the 
modified  spheroid  are  given  in  Table  2. 


Source 

CD 

Experiment  - Patel  et  al.  [15] 

Present  Calculation 

Calculation  - Nakayama  et  al.  [10] 

Calculation  - Myring  [6] 

.00447  (average) 

.00454 

.00485 

.00464 

Table  2.  Drag  Coefficient  for  Modified  Spheroid. 
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The  present  method  gives  the  value  closest  to  experiment. 

Note  that  Patel  and  Guven  [12]  show  that  the  use  of  Granville's 
wake  formula,  Eq.  (24),  leads  to  an  overprediction  of  for  the 

spheroid  of  30  percent  where  presumably  they  are  using  measured  values 
of  0 and  at  the  tail  station.  The  present  prediction  of  CD  does 

not  bear  out  their  finding  possibly  because  Weinbaum's  hypothesis  is  used 
here  to  determine  u^  at  the  tail  and  0 is  determined  differently. 

F-57  Body 

The  low  drag  F-57  body  tested  by  Patel  and  Lee  [15]  is  the  last  case 
considered.  This  body,  shown  in  Fig.  16,  has  a nearly  cusped  tail  so  that 
the  boundary  layer  on  the  aft  end  has  a different  history  from  the  previous 
cases.  The  chord  Reynolds  number  of  the  experiment  was  1.20  x 10^  and  the 
boundary  layer  was  tripped  at  x^  = 0.475,  just  ahead  of  where  laminar 
separation  would  have  occurred.  Comparisons  between  present  calculations, 
those  of  Patel  and  Lee  [15]  and  experiment  are  shown  in  Figs.  17  - 23. 

The  computed  results  of  Patel  and  Lee,  which  are  based  on  the  integral 
boundary-layer  method  of  Patel  [17],  used  the  measured  pressure  distribution 
at  the  edge  of  the  boundary  layer  plus  estimated  values  of  the  pressure 
integrals  arising  from  variations  normal  to  the  streamlines.  Integration 
of  the  equations  was  started  at  x^  = 0.70.  Hence  the  calculations  of 
Patel  and  Lee  are  not  a true  strong  interaction  solution  in  which  the 
pressure  distribution  is  part  of  the  solution. 

As  Fig.  17  shows,  the  present  pressure  distribution  calculation  agrees 
well  with  experiment  both  on  the  body  and  in  the  wake.  As  in  previous  cases 
the  calculated  result  is  slightly  higher  than  experiment  in  the  peak  region 
near  the  tail.  For  this  body  comparison  with  the  potential  flow  pressure 
distribution  shows  that  the  strong  interaction  region  begins  at  = 0.75. 


r 
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In  computing  the  boundary-layer  solution  for  the  F-57  body,  the  station 
where’  the  eddy  viscosity  is  introduced  was  varied  until  agreement  was  obtained 
with  the  measured  momentum  deficit  area  in  the  trip  region.  As  Fig.  18  shows, 
the  calculated  results  for  0 by  both  methods  agree  well  with  experiment  on  the 
body  but  agreement  becomes  increasingly  worse  downstream  in  the  wake.  Conse- 
quently a question  arises  as  to  the  source  of  the  discrepancy  between  calculation 
and  experiment  for  0 in  the  wake.  Estimates  by  Patel  and  Lee  (see  their  Fig.  26) 
indicate  that  the  pressure  integral  Ip  arising  from  streamline  curvature  is 
practically  negligible  beyond  the  near  wake  (Xq  greater  than  about  1.04). 

Hence  the  wake  momentum  integral  equation,  Eq.  (19),  should  adequately  represent 
the  variation  of  0,  especially  if  the  measured  values  of  H and  u^  (corresponding 

to  C on  the  wake  centerline)  are  used.  Such  a calculation  was  made  by  numeri- 
P 

cally  integrating  Eq.  (19)  beginning  from  measured  values  at  x^  = 1.04.  The 
results,  shown  in  Table  3,  are  practically  the  same  as  the  present  strong  inter- 
action solution  indicating  that  the  experimental  values  of  0 in  the  wake  by 
Patel  and  Lee  are  in  error.  Consequently,  the  drag  coefficient  given  by  Patel 


x0/L 

U 

e 

exper . 

H 

exper . 

0 x 104 
exper. 

0 x 104 
calc . (1) 

0 x 104 
calc . (2) 

1.04 

0.9690 

1.319 

1.509 

1.509 

1.521 

1.06 

0.9726 

1.305 

1.464 

1.491 

1.497 

1.10 

0.9803 

1.279 

1.361 

1.453 

1.446 

1.20 

0.9899 

1.262 

1.291 

1.407 

1.402 

1.30 

0.9915 

1.235 

1.197 

1.400 

1.390 

1.40 

0.9965 

1.222 

1.158 

1.377 

1.375 

2.47 

0.9999 

1.124 

1.007 

1.362 

CO 

1.0000 

1.000 

1.361 

1.354 

1 - Momentum  eauation  with  exoerimental  H and  u . 

e 

2 - Present  strong  interaction  method. 

Table  3.  Variation  of  Wake  Integral  Parameters  for  F-57  Body. 
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and  Lee  (based  on  maximum  frontal  area)  is  too  low.  In  fact  their  value, 
based  on  the  measured  0 at  x^  = 2.47,  is  incorrect.  They  give  Cg  = 0.0092 
whereas  it  should  be  0.0294.  The  present  calculation  gives  CQ  = 0.0396. 

The  large  trip  wire  used  in  the  experiment  (1.664  mm  in  diameter)  is 
responsible  for  the  measured  values  of  C^,  shown  in  Fig.  19,  being  considerably 
higher  than  the  computed  curve  immediately  downstream  of  the  transition  region. 
Beyond  about  x^  = 0.70  the  effect  of  the  wire  dies  out  and  the  two  calculated 
curves  agree  well  with  experiment.  The  trip  wire  also  manifests  itself  in 
the  shape  factor  distribution  as  shown  in  Fig.  20.  Over  the  last  20  percent 
of  body  length  the  present  calculation  is  too  low  but  agrees  well  with 
experiment  in  the  wake.  The  cause  of  the  low  H prediction  on  the  aft 
end  of  the  body  is  the  algebraic  turbulence  model,  as  in  the  previous  two 
cases.  The  shape  factor  is  predicted  by  the  present  method  more  accurately 
for  a concave  tail  than  a convex  one.  The  velocity  profiles,  shown  in 
Figs.  21  - 23,  illustrate  respectively  the  trip  wire  effect  near  the  wall 
at  Xq  = 0.601,  the  usual  overprediction  of  the  velocity  near  the  wall  at 
Xq  = 0.88  and  a spurious  bump  near  the  wall  in  the  calculated  profile  at 
Xq  = 0.99  caused  by  extreme  stretching  of  the  normal  step  size  from  the 
Mangier  transformation  singularity  at  the  tail. 
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6 . Conclusions 

The  pressure  mapping  hypothesis  of  Weinbaum  has  been  found  to  provide 
a simple  means  of  adapting  the  displacement  body  concept  to  simulate  stream- 
line curvature  effects  in  the  high  Reynolds  number  axisymmetric  strong 
interaction  problem.  Calculated  results  by  the  present  method  for  the  three 
cases  considered  generally  agree  well  with  experiment  except  for  the  boundary- 
layer  shape  factor.  The  predictions  by  the  present  method,  where  a finite 
difference  boundary-layer  solution  is  used,  are  about  the  same  in  accuracy 
as  the  predictions  of  a good  integral  boundary-layer  procedure  such  as  that 
of  Patel  [17]  , except  for  the  boundary-layer  shape  factor  which  is  predicted 
better  by  the  integral  method.  The  reason  for  the  poor  prediction  of  H 
in  the  tail  region  by  the  present  method  is  the  basic  inadequacy  of  an 
albebraic  turbulence  model  in  that  region.  Such  a model  does  not  properly 
represent  the  change  from  a boundary-layer  profile  to  a wake  profile  nor 
does  it  include  streamline  curvature  effects.  The  inadequacy  of  algebraic 
turbulence  models  in  this  regard  has  been  noted  by  other  investigators 
[3,  14,  15,  18]. 

The  present  calculated  results  for  body  drag  coefficient  have  been  found 
to  agree  well  with  experiment  and  to  be  as  good  or  better  as  other  published 
calculated  values,  at  least  for  the  Lyon  model  A body  and  the  modified 
spheroid  tested  by  Patel  et  al.  For  the  case  of  the  low  drag  F-57  body  some 
doubt  exists  as  to  the  validity  of  the  measured  momentum  deficit  area  in  the 
wake.  The  simple  integral  method  of  Granville  has  been  used  to  model  the  wake 
development  and  hence  predict  0 which  is  sensitive  to  the  initial  value  of 

oo 

u^  in  the  wake.  By  using  Weinbaum' s hypothesis  u^  has  been  corrected  to 
account  for  streamline  curvature  effects  in  the  tail-near  wake  region  so 

that  the  predicted  value  of  0 , and  hence  C , is  reasonably  accurate. 

OO  D 
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.Improvements  in  the  present  calculation  method  could  be  made  in  the 

following  areas,  in  addition  to  modifying  the  turbulence  model: 

1.  The  approximation  of  concentric  streamlines,  Eq.  (9),  which  leads  to 
the  simple  form  of  Weinbaum's  hypothesis,  could  be  changed  to  a more 
realistic  assumption.  The  resulting  expression  for  the  pressure  change 
across  the  boundary  layer  would  undoubtedly  lead  to  a pressure  mapping 
more  complicated  than  that  of  Weinbaum. 

2.  The  boundary-layer  equations  could  be  integrated  in  finite  difference 
form  downstream  in  the  wake,  as  well  as  along  the  body,  to  eliminate 
the  incompatibility  of  a finite  difference  solution  with  an  integral 
solution.  As  long  as  boundary-layer  coordinates  are  retained,  this 
approach  would  work  properly  only  for  a body  with  a cusped  tail.  For 
a body  with  a nonzero  tail  angle  boundary-layer  coordinates  do  not 
mesh  properly  with  cylindrical  wake  coordinates  resulting  in  an 
overlap  region  in  the  coordinate  systems.  To  restart  the  numerical 
calculation  in  the  wake  would  require  interpolation  of  the  boundary- 
layer  solution.  Such  an  interpolation  would  not  be  satisfactory 
because  the  boundary-layer  approximation  is  coordinate  system 
dependent.  In  addition,  the  Cebici-Smith  boundary-layer  scheme 
cannot  be  used  in  the  wake  because  of  the  singularity  at  zero  radius 
in  the  Mangier  transformation.  Physical  coordinates  should  be  used 
for  a finite  difference  wake  calculation. 
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One  possibility  of  overcoming  the  coordinate  mismatch  problem  at  the 
tail-wake  juncture  is  to  abandon  the  usual  boundary-layer  coordinate  system 
and  use  either  a conformal  system  or  a skewed  system.  At  this  stage  the 
artificial  division  of  the  flow  into  a boundary  layer  and  an  outer  inviscid 
region  hardly  seems  worthwhile.  In  ".he  new  coordinates  just  mentioned  one 
could  proceed  to  a parabolized  form  of  the  Navier-Stokes  equations,  with 
an  appropriate  turbulence  model,  to  be  solved  iteratively  for  the  entire 
strong  interaction  region  similar  to  the  scheme  of  Schetz  and  Favin  [2] . 

If  done  properly  the  machine  time  required  to  solve  such  a system  should 
be  little  great(A-  than  the  strong  interaction  procedure  described  herein. 

Such  an  alternative  is  presently  under  investigation. 
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Figure  1.  Displacement  Body  Geometry 


Figure  2.  Lyon  Model  A Body  Geometry 
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Figure  4.  Momentum  Deficit  Area,  Lyon  Model  A Body 
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157.8  cm 

BOUNDARY-LAYER  TRIP 


Figure  9.  Modified  Spheroid  Geometry 
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Figure  11.  Momentum  Deficit  Area,  Modified  Spheroid 
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Figure  20.  Boundary-Layer  Shape  Factor,  F-57  Body 
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